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Chromospheric Polarization 
in the Photospheric Solar Oxygen Infrared Triplet 

Tanausu Del Pino Aleman^’^ & Javier Trujillo Bueno^’^’^ 

ABSTRACT 

We present multilevel radiative transfer modeling of the scattering polarization 
observed in the solar O i infrared triplet around 777 nm. We demonstrate that the 
scattering polarization pattern observed on the solar disk forms in the chromosphere, 
far above the photospheric region where the bulk of the emergent intensity profiles 
originates. We study the sensitivity of the polarization pattern to the thermal structure 
of the solar atmosphere and to the presence of weak magnetic fields (10“^ - 100 G) 
through the Hanle effect, showing that the scattering polarization signals of the oxygen 
infrared triplet encode information on the magnetism of the solar chromosphere. 

Subject headings: line: profiles — polarization — scattering — radiative transfer — 
Sun: chromosphere — Sun: photosphere — Stars: atmospheres 


1. Introduction 


The absorption and scattering of anisotropic radiation in the solar atmosphere produce a richly 
structured lin early polarized spectrum, which is of great diagnostic potential (e.g., the review by 
Stenfioll2009ll . In general, the physical origin of this spectrum is the radiatively induced population 
imbalances and quantum coherence between the magnetic sublevels of the line’s levels (i.e., atomic 
level polarization). A weak magnetic field (i.e., such that the level’s Zeeman splitting equals the 
inverse of its lifetime) is sufficient to produce a significant modification of the atomic level polar¬ 
ization (Hanle effect). In contrast to the Hanle effect, elastic collisions may destroy completely 
the atomic level polarization if their rates are significantly larger than the inverse of the level’s 
lifetime. Since ground and metastable levels are relatively long-lived, they are more vulnerable 
than excited levels to the above-mentio ned effects of magnetic fields and elastic collisions (e.g., 
Landi Degl’Innocenti & Landolfi 2004). 
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The scattering polarization observed in a speetral line has, in general, two eontributions: one 
due to the seleetive emission of polarization eomponents, eaused by the presenee of atomie po¬ 
larization in the upper level, and another due to seleetive abso rption of polarization eomponents. 


eaused by the presenee of atomie polarization in the lower level ( 

Truiillo Bueno & Landi Deal’Innoeenti 

1997 

Truiillo Buendl999: 

Truiillo Bueno et al.ll2002: Manso Sainz & Truiillo BuendboOSa 

). That 


the atomie polarization of the lower level ean eontribute to the observed seattering line polariza¬ 
tion ean be easily understood by eonsidering a one-dimensional (ID) model atmosphere, whose 
plasma ean be unmagnetized or magnetized by a mieroturbulent and isotropieally distributed weak 
magnetie field. In both oases the transfer equation for the Stokes parameter Q reads (with s the 
geometrioal distanoe along the ray): 


ds 


[eg - TJqI] - T]iQ, 


( 1 ) 


where eg and rjx (with X = I or X = Q) are the eorresponding emission and absorption ooeflioients. 
The absorption ooeffioient rjg is proportional to the alignment Pq(^) of the lower level (whioh quan¬ 
tifies the population imbalanoes among the lower level sublevels), while the emission eoeffieient 
eg is proportional to the alignment of the upper level pliu) (whieh quantifies the population imbal- 
anees between the sublevels of the upper level). Sinee in Eq. (1) the souree of Stokes Q is given 
by the term within braekets, we have that Stokes Q ean be produeed either beeause the emission 
eoeffieient €q^0 (seleetive emission of polarization eomponents, whieh requires and/or 

beeause the absorption eoeffieient rigi^O (seleetive absorption of polarization eomponents, whieh 
requires Pq(/’)?^ 0). This last eo ntribution is termed “zero-field” diehro ism, beeause it has nothing to 
do with the Zeeman effeet (ef.. lManso Sainz & Trujillo Buenoil2003ah . If elastie eollisions are very 
effieient 6g «? 7 g « 0, while if they are moderately effieient we may have and no diehroism at 
all (i.e., riQ^O). In general, both eg and rjg may be signifieant. 

An example of solar speetral lines whose observed fraetional linear polarization Qjl profiles 
show elear signatures of sueh t wo eontributions is that of the IR triple t of Ca ii observed in quiet re¬ 
gions of the solar atmosphere (IManso Sainz & Trujillo Buenoll2003ah . Here elastie eollisions with 
neutral hydrogen atoms do not destroy the polarization of the line’s upper and lower levels in the 
ehromospherie region where the eore of the Stokes I profiles originates. Under sueh eireumstanees 
the following approximate expression ean be used to obtain an estimatffi n of the amplitude of the 
fraetional seattering polarization at the line eenter (ITrujillo Buenoll2003h 
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crlUt)], 


( 2 ) 


where w®, is a eoeffieient introdueed bv iLandi Degl’Innoeentil (119841) whose value depends on the 


total angular momentum of the line’s upper and lower levels (see Table 10.1 of lUandi Degl’Innoeenti & Landolfi 
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Fig. 1.— Examples of the observations by iTrujillo Bueno et al.l (1200111 of the seattering polariza¬ 
tion in the IR triplet of O i at the three indieated on-disk positions from the nearest limb. The 
positive referenee direetion for Stokes Q is the parallel to the nearest limb. 


20041) and cr\{J) = (with Po(7) proportional to the overall population of the 7-level) 

quantifies the level’s fraetional atomie alignment at the height in the model atmosphere where the 
line eenter optieal depth is unity along the line of sight (LOS). The 0-^(7) values have to be obtained 
by solving the statistieal equilibrium equations for the multipole moments of the atomie density 
matrix assoeiated to eaeh atomie level 7 and the Stokes-veetor transfer equation for eaeh of the 
radiative transitions in the atomie model. 

Speetropolarimetrie observations of the quiet solar atmosphe re in the infrared (IR) tripl et of 
O I show ele ^ seattering polarizatio ns signals, both on the disk (ITrujillo Bueno et al.ll200lh and 
o ff the limb dSheelev & Keller! 120031) . Figure [H adapted from a similar one in the review paper 
of 


me timb (laneetev & KettenizUUil) . PigureLLl adapted trom a simitar one m me review paper 
Trujillo Buenol (200^, s hows some elose to the limb examples of the on-disk observations by 


Trujillo Bueno et al.l (120011) . Note that while the lines at 7772 A (hereafter, line 1) and 7774 A 
(hereafter, line 2) showed positive Q/l signals, with line 2 slightly more polarized than line 1, 
negative Qjl values were fou nd for the 7775 A line (hereafter, line 3). Given the eomplexity of 
the radiative transfer problem. IBelluzzi & Trujillo Buenol (1201 ih assumed a plane-parallel layer of 
O I atoms illuminated by the solar eontinuum radiation field, and applied Eg . (2) a fter solving the 


statistieal equilibrium equations for the values (see also lTrujillo Buenoll2009f) . These authors 


eoneluded that the observed Qjl signals ean be qualitatively explained if there is a sizable amount 
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of atomic polarization in the upper and lower line levels around the height in the solar photosphere 
where the optieal depth at the line eenter is unity along the line of sight, with zero-field diehroism 
playing a signifieant role. 

However, as shown in this letter, the situation is more subtle and interesting. It turns out that in 
the solar photosphere elastic collisions between the O i and H i atoms destroy the atomie level po¬ 
larization of the O IIR triplet (i.e., crliJe) » (t^{Ju) ~ 0 in the photospheric region where the bulk of 
the emergent intensity profiles originates, both for the lower and upper levels). Therefore, aecord- 
ing to the approximate Eq. (2) the O i IR triplet should show negligible polarization, eontrary to the 
on-disk observations of Fig. 1. In order to find the true physieal origin of the observed seattering 
polarization, we solved numerieally the full (non-loeal and non-linear) radiative transfer problem 
iteratively and ealeulated the emergent Stokes I and Q profiles using the resulting self-eonsistent 
values of the emission and absorption eoeffieients at eaeh point along the line of sight. The main 
result we want to highlight here is that the seattering polarization signals observed in the oxygen 
IR triplet are produced in the solar chromosphere. As we shall see below, their physieal origin is 
the seleetive emission and selective absorption (diehroism) processes eaused by the presenee of 
atomic level polarization within the solar chromosphere. 


2. Theoretical approach 


Our theoretieal approach is based on th e quantum-meehanieal density matrix theory explained 
in Landi DeglTnnoeenti & Landolfil (l2004h . which assumes that the absorption and re-emission of 
photons are statistically independent events. Belluzzi & Trujillo Bueno ( 201ll) have investigated 
the impact of quantum interference between pairs of sublevels pertaining to different fine-structure 
y-levels (hereafter, 7-state interference), showing that it can be neglected for modeling on-disk 
observations of the scattering polarization in the O i IR triplet. For the case of a multilevel atom 
neglecting 7-state interference the above-mentioned theory is a reliable approach whenever the 
radiation field that pumps the atomic system is flat over a frequency interval larger than the nat¬ 
ural width of the levels. This is a suitable assumption for modeling the oxygen IR triplet and 
we have developed a radiative transfer code that solves jointly the statistical equilibrium equa¬ 
tions for the multipole moments of the atomic density matrix associated to each atomic level 7 
and the Stokes-vector transfer equation for each of the radiative transitions in the atomic model . 
The method of solution is based on a Jacobian multilevel iterative scheme ( Trujillo Buenol 1999 : 
Manso Sainz & Trujillo Buenql i2003bh and on the BESSER formal solver of the Stokes-vector 
transfer equation proposed bv IStepan & Trujillo Buenol (l2013h . which we have used also for com¬ 
puting the emergent Stokes I and Q profiles that give the Q/I ones shown below. 


The atomic model, which has 23 O i 7-levels and the ground level of O ii, is an improved ver- 
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sion of the model atom used bv ICarlsson & Judge! (119931) . The model takes into acc ount 43 bound- 


bound radiative transitions, with Einstein coefficients taken from the NIST database (iKramida et al. 


2013h . The number density of neutral oxygen atoms at each atmospheric height was calculated 


taking into account the ionizations from each O i level, with the cross sections taken from the 
Opacity ProjecR The bound-bound inelastic collisions are calculated using the approximation of 


van Regemorten (Il962h if the corre s pondi ng transition is radiatively allowed and with the approx- 

V forbidden . 

M (l200nl). 


Bely & van Regemorten ( 


19701) if the corresponding transition is radiative 


imation of 

The bound-free inelastic collisions are calculated following the formulae given by 
Collis ions with protons are taken into account for the ground term levels, following 
(Il977l). 


Haisch et al. 


The lines of the O i IR triplet share the same (metastable) lower level, whose angular momen¬ 
tum is 7^ = 2; the angular momentum values of the three upper levels of increasing energy are 1, 
2 and 3. Therefore all these levels can carry atomic alignment and contribute to the observed scat¬ 
tering polarization. When observed on the solar disk at a line of sight with /r = cos6 = 0.1 (6 being 
the heliocentric angle), the mean formation heights of the line core intensity profiles are located 
around 200 km above the solar surface. At such low photospheric heights the number density of 
neutral hydrogen atoms is of the order of 10'^cm~^, so we have to investigate if elastic collisions 
between neutral oxygen and hydrogen atoms are able to depolarize the Q i IR tri plet levels. To this 
end, we have used equation (7.108) of iLandi DeglTnnocenti & Landolfil (120041) . which is suitable 
to estimate the rates D® of such elastic collisions for the lower (z = £) and upper (z = u) levels of 
the O I IR triplet. 


Figure [2] shows the variation with height in a semi-empirical model of the quiet solar atmo¬ 
sphere of the rates of elastic collisions with neutral hydrogen, for the lower level and the three 
upper levels of the O i triplet. These collisional rates have to be compared with the rates of spon¬ 
taneous radiative emissions (upper dashed line) and with the rates of radiative absorptions (lower 
dashed curve). As can be seen, below 400 km both radiative rates are smaller than the rates of 
elastic collisions and, therefore, we may expect that in the low photosphere of the quiet Sun (see 
the shaded region in Fig. ^ elastic collisions with neutral hydrogen atoms are efficient in destroy¬ 
ing the atomic alignment that anisotropic radiation pumping processes induce in the O i IR triplet 
levels. 


'http://vizier.u- strasbg.fr/topbase/topbase.html 
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Fig. 2.— The variation with height in a ID semi-empirieal model of the quiet solar atmosphere 
of the rates of elastie eollisions with H i atoms, for the levels of the O i IR triplet (solid eurves). 
The upper and lower dashed lines indieate the height variation of the rates of spontaneous ra¬ 
diative emissions and of radiative absorptions, respeetively. The shaded area indieates the atmo- 
spherie heights where the optieal depth at the eenter of the O i lines is unity for line of sights with 


3. Results 

This seetion shows and diseusses the results we have obtained by solving jointly the statistieal 
equilibrium and radiative transfer equations, using the above-mentioned multilevel oxygen model 
and semi-empirieal models of the quiet solar atmosphere. 


3.1. The atomic polarization and the emergent spectral line polarization 


Figure |3] shows the height variation of the self-eonsistent fraetional atomie alignment values 
of the O I IR triplet levels. The solid eurves take into aeeount the impaet of elastie eollisions with 
neutral hydrogen atoms, while the dashed eurves eorrespond to the ease in whieh elastie eollisions 
are negleeted. To obtain sueh = pl(J)lp^(J) values we numerieally solv ed the ensuing 


radiative transfer problem of seattering line polarization in the quiet Sun model C of lFontenla et al. 


(119931) (hereafter, FAL-C model) using the above-mentioned realistie multilevel oxygen model and 
taking into aeeount the impaet of the hydroge n Lyman-/? pumping on the populations of the oxygen 
levels (see iMiller-Rieei & Uitenbroek 2002 ). As shown by the solid eurves of the figure, in the 
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model’s photosphere elastie eollisions strongly reduee the polarization of the O i IR triplet levels, 
yielding (t^{J) « 0 in the photospherie formation region of the Stokes-/ profiles. Therefore, sinee 
this implies that 6 q« 0 and ?7 q« 0 in the solar photosphere, the observed Q/I signals shown in Fig. [T] 
eannot be produeed in the same atmospherie region (the photosphere) where the bulk of the Stokes 
I profiles originates. 

Interestingly, if elastie eollisions are fully neglected (see the dashed curves of Fig. [3]) the 
ensuing Q/l profiles of the calculated emergent radiation do not agree with the observations of 
Fig. [T](see the dashed curves of the central panels of Fig. [Hand note that without elastic collisions 
the calculated Q!I signal for line 3 is positive, while the observed one is negative). However, as 
shown by the solid curves of the central panels of Fig. HI the calculated Q!I profiles qualitatively 
agree with the observed ones as soon as the impact of the elastic collisions is taken into account. 


In retrospect, we can clarify the reason why the application of the app r oxima te Eddington- 
Barbier Eq. (2) used by Trujillo Bueno ( 2009 ) and Belluzzi & Trujillo Bueno ( 201 ih allowed them 
to explain qualitatively the observations of Eig. 1, although missing the interesting point explained 
in this Eetter (i.e., that the observed scattering polarization originates in the chromosphere). As 
mentioned in Section 1, those authors considered a plane-parallel layer of O i atoms and solved 
the statistical equilibrium equations for obtaining the cr^iJ) values assuming that the anisotropic 
pumping is exclusively due to the continuum radiation coming from the underlying solar visible 
surface. Such a predominantly vertical incident radiation is characterized by a positive anisotropy 
factor, which in the absence of elastic collisions leads to significant crp(7) values with signs similar 
to those found in Eig. 3 at chromospheric heights. 


3.2. Zero-field chromospheric dichroism 

In order to demonstrate that zero-field dichroism within the solar chromosphere produces clear 
observational signatures, we have calculated the 2/7 profiles that result when the contribution of 
the tjqI term of Eq. (1) is neglected above 500 km in the EAE-C semi-empirical model. The dashed 
curves of the bottom panels of Eig. Hshow the ensuing results, along with the previously shown 
2/7 profiles (solid curves of Eig. H]) that take into account the full contribution. Clearly, the Q/I 
profile of line 3 is dominated by zero-field chromospheric dichroism. The reason why the Q/I 
signal of line 3 is dominated by the tiqI term of Eq. (1) (i.e., by zero-field dichroism) is easy to 
understand by noting that for line 3 (see Eig. [3]), while wf^j ~0.6 and = 0.1. 

Moreover, zero-field dichroism has a significant impact on the fractional polarization line-center 
amplitudes of lines 1 and 2. As seen in the bottom panels of Eig. HI only when the full impact of 
zero-field dichroism is taken into account we obtain that line 1 is slightly less polarized than line 
2, in agreement with the observations of Eig. [T] 
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3.3. Sensitivity to the model’s thermal structure 


The scattering polarization depends on the anisotropy of the radiation field that produces 
optical pumping within the solar atmosphere, which in turn depends on the thermal structure of the 
atmospheric plasma among other factors. For this reason, we have investigated also the sensitivity 
of the emergent Q/l profiles of the O i IR triplet to the thermal structure of the model atmosphere. 
To that end, we carried out radiative transfer calc ulations in tw o semi-empirical models: the FAL-C 
model already considered and the Mco m odel of|Avretti(ll995f) which is signific antly cooler around 
a height of 1000 km (e.g., see figure [2l of IManso Sainz & Trujillo Buenoll2010h . 


Since the two atmospheric models are practically identical below 500 km, we point out that 
the emergent Stokes-/ profiles are similar. However, the calculated 2// profiles turn out to be very 
different, in agreement with the previously reported point that the scattering polarization signals are 
produced within the chromosphere where the two models clearly differ. Interestingly enough, only 
the Qjl profiles calculated in the (hotter) FAL-C model agree qualitatively with the observations 
of Fig. [I](e.g., we find that in the Mco model line 3 has a positive signal, contrary to the observed 
Q!I profile). This additional result emphasizes that the observed scattering polarization is also of 
interest for constraining the thermal structure of the solar atmosphere. 


3.4. Hanle sensitivity to the presence of a magnetic field 

The critical magnetic field Rh (m gauss) needed to produce a significant modification of the 
polarization of an atomic level via the Hanle effect is that for which the level’s Zeeman splitting (in 
frequency units) equals the inverse of its radiative lifetime: fin = (l-137xl0“^)/(gtiife) (with g the 
level’s Lande factor and tufe its radiative lifetime in seconds). Since the lower level of the O i IR 
triplet is metastable (i.e., long-lived), milligauss fields are sufficient to significantly alter its atomic 
polarization, while the critical fields of the three upper levels of the O i IR triplet are 1.7 G, 2.2 G 
and 2.5 G for the upper levels of lines 3, 2 and 1, respectively. 

Figure [5] illustrates the sensitivity to the Hanle effect of the scattering polarization signal at the 
center of the O i IR triplet lines, assuming random-azimuth magnetic fields of fixed inclination. As 
shown in the figure, the scattering polarization of these O i lines react to magnetic fields between 
10“^ and 100 gauss. Finally, we point out that the response functions of the line-center amplitudes 
of the Q/l profiles to magnetic strength perturbations in the FAL-C semi-empirical model show 
significant values between 600 and 1600 km, approximately, with their peaks around 1000 km . 
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4. Concluding comments 

The intensity profiles of the solar O i IR triplet originate mainly in the low photosphere. How¬ 
ever, we have shown that their seattering polarization signals are generated above the photosphere, 
where elastie eollisions with neutral hydrogen atoms do not destroy the atomie polarization that 
anisotropie radiation pumping proeesses induee in the oxygen levels. The physieal origin of the 
observed linear polarization signals is the seleetive absorption (diehroism) and seleetive emission 
of polarization eomponents eaused by the atomie level polarization of the ehromosphere. The 
seattering polarization of the O i IR triplet is sensitive to the thermal and magnetie strueture of 
the ehromosphere. In partieular, via the Hanle effeet the seattering polarization amplitudes are 
sensitive to magnetie fields with strengths between milligauss and a few gauss, and with response 
funetions to magnetie field perturbations that peak around 1000 km above the model’s visible sur- 
faee. 


New speetropolarimetrie observations of the O i IR triplet would be very weleomed. Of 
speeial interest would be to observe regions with different levels of magnetie aetivity (e.g., within 
and outside eoronal holes) and to pay attention to the sign of the seattering polarization signals in 
the three lines and to their relative polarization amplitudes. By eontrasting sueh observables with 
radiative transfer ealeulations in inereasingly realistie models of the solar atmosphere we may hope 
to learn more about the thermal and magnetie strueture of its enigmatie ehromosphere. 


We are grateful to Rafael Manso Sainz (lAC) for helpful diseussions and suggestions. Finan- 
eial support by the Spanish Ministry of Eeonomy and Competitiveness through projeets AYA2010- 
18029, AYA2014-60476-P and CONSOLIDERINGENIO CSD2009-00038 are gratefully aeknowl- 
edged. 
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Fig. 3.— The height variation of the fraetional atomie alignment of the O i IR triplet levels ealeu- 
lated in the FAL-C semi-empirieal model. The solid eurves take into aeeount the elastie eollisions 
with neutral hydrogen atoms while the dashed eurves negleet them. Note that the top horizontal 
axes give the line eenter optieal depth of line 1 along a LOS with jj. = 0.1. 
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Fig. 4.— The Stokes 1 and Q/l profiles of the IR triplet of O i calculated in the FAL-C semi- 
empirical model for a line of sight with jj. = 0.1. The positive reference direction for Stokes Q 
is the parallel to the nearest limb. Top panels: the Stokes I profiles normalized to the continuum 
intensity. Central panels: Q/I profiles calculated with (solid curves) and without (dashed curves) 
elastic collisions. Bottom panels: Qjl profiles neglecting the contribution of dichroism (dashed 
curves) and taking it into account (solid curves). 
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Line 1 Line 2 Line 3 



Fig. 5.— The Qjl line-center amplitudes of the IR triplet of O i versus the strength of a micro- 
structured magnetic field with random azimuth and the indicated inclination with respect to the 
local vertical, calculated in the FAL-C semi-empirical model for a line of sight with /j. = 0.1. The 
positive reference direction for Stokes Q is the parallel to the nearest limb. 





